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Abstract 

In this work, we consider rational ordinary differential equations dy/dx = 
Q(x, y)/P(x, y), with Q(x, y) and P(x, y) coprime polynomials with real coef- 
ficients. We give a method to construct equations of this type for which a first 
integral can be expressed from two independent solutions of a second-order 
homogeneous linear differential equation. This first integral is, in general, 
given by a non Liouvillian function. 

We show that all the known families of quadratic systems with an irre- 
ducible invariant algebraic curve of arbitrarily high degree and without a ra- 
tional first integral can be constructed by using this method. We also present 
a new example of this kind of families. 

We give an analogous method for constructing rational equations but by 
means of a linear differential equation of first order. 
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1 Introduction 



This paper deals with rational ordinary differential equations such as 

dy_ Q(x,y) 

dx P(x,y) ' 1 ' 

where Q(x,y) and P(x,y) are coprime polynomials with real coefficients. We asso- 
ciate to this rational equation a planar polynomial differential system by introducing 
an independent variable t usually called time. Denoting by ' = d/dt, we have 

x = P(x,y), y = Q(x,y), (2) 

where {x,y) G 1R 2 . This system defines the vector field X = P{x,y)^ + Q(x,y)j^ 
over M? and, equivalently, the 1-formo; = Q(x,y)dx — P(x,y)dy. We indistinctively 
talk about equation (JTJ) and system (0). Let d be the maximum degree of P and Q. 
We say that system © is of degree d. When d = 2, we say that © is a quadratic 
system. 

In order to simplify notation, we define M[x, y] as the ring of polynomials in 
two variables with real coefficients and lR(x, y) as the field of rational functions in 
two variables with real coefficients, that is, the quotient field of the previous ring. 
Analogous definitions stand for R[x) and WL(x). 

We have an equation (0) defined in a certain class of functions, in this case, the 
rational functions with real coefficients M.(x, y) and we consider the problem whether 
there is a first integral in another, possibly larger, class. For instance, as we will 
discuss later on, H. Poincare stated the problem of determining when a system (J2J) 
has a rational first integral. A C k function H : U — > K such that it is constant 
on each trajectory of (0) and it is not locally constant is called a first integral of 
system (J2]) of class k and the equation H(x,y) = c for a fixed c6l gives a set of 
trajectories of the system, but in an implicit way. When k > 1, these conditions 
are equivalent to w A dH = and H not locally constant. The problem of finding 
such a first integral and the functional class it must belong to is what we call the 
integrability problem. 

To find an integrating factor or an inverse integrating factor for system (0) is 
closely related to finding a first integral for it. When considering the integrability 
problem we are also addressed to study whether an (inverse) integrating factor 
belongs to a certain given class of functions. 

Definition 1 Let W be an open set ofM 2 . A function fi : W — > R of class C k (W), 
k > 1, that satisfies the linear partial differential equation 

to A dfi = fi du, (3) 

is called an integrating factor of system fflj) on W. 
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It has been shown than an easier function to find which also gives additional prop- 
erties for a differential system J2J is the inverse of an integrating factor, that is, 
V — which is called inverse integrating factor. 

We note that {V = 0} is formed by orbits of system (j2J). The function fi = 1/V 
defines on W \ {V = 0} an integrating factor of system (j2J), which allows the 
computation of a first integral of the system on W \ {V = 0}. The first integral H 
associated to the inverse integrating factor V can be computed through the integral 
H(x,y) = J u/V, and the condition Q for /i = 1/V ensures that this line integral 
is well defined. 

The inverse integrating factors play an important role in two of the most difficult 
open problems of qualitative theory of planar polynomial vector fields, which are 
the center problem and 16 th Hilbert problem. In jB], it has been noticed that for 
many polynomial differential systems with a center at the origin there is always an 
inverse integrating factor V globally defined in all M. 2 , which is usually a polynomial. 
However, the first integral for a polynomial differential system with a center at the 
origin can be very complicated. 

We say that a function f(x, y) is an invariant for a system (j2J) if uAdf = kf with 
k(x,y) a polynomial of degree lower or equal than d — 1, where d is the degree of 
the system. This polynomial k(x, y) is called the cof actor oi f(x, y). In the previous 
equality and all along this paper we use the convention of identifying the space of 
functions over M 2 and the space of 2-forms over M. 2 . In case f(x,y) = defines a 
curve in the real plane, this definition implies that u A df equals zero on the points 
such that f(x,y) = 0. In case f(x,y) is a polynomial we say that f(x,y) = is an 
invariant algebraic curve for system (j2J). 

Let us consider f(x,y) = an invariant algebraic curve for system (j2J), we will 
always assume that f(x, y) is an irreducible polynomial in M.[x, y]. Otherwise, it can 
be shown that each of its factors is an invariant algebraic curve for system (j2J). We 
will denote by n the degree of the polynomial f(x,y). 

In [T2] , G. Darboux gives a method for finding an explicit first integral for a 
system (0) in case that d(d + l)/2 + l different irreducible invariant algebraic curves 
are known, where d is the degree of the system. In this first integral of 

the form H = f^f^ 2 ■ ■ ■ fr r , where fi(x,y) = is an invariant algebraic curve for 
system (j2J) and A,, G C not all of them null, for i = 1, 2, . . . , r, r 6 N, can be defined 
in the open set R 2 \ S, where S = {(x, y) E R 2 | (A • f 2 ■ . . . • f r )(x, y) = 0}. The 
functions of this type are called Darboux functions. We remark that, particularly, if 
\i G Z , Vi = 1, 2, . . . , r, H is a rational first integral for system (j2J). In this sense J. 
P. Jouanoulou [TH], showed that if at least d(d+ 1) +2 different irreducible invariant 
algebraic curves are known, then there exists a rational first integral. 

The main fact used to prove Darboux's theorem (and Jouanoulou's improve- 
ment) is that the cofactor corresponding to each invariant algebraic curve is a 
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polynomial of degree < d — 1. Invariant functions can also be used in order to 
find a first integral for the system. This observation permits a generalization of 
Darboux's theory which is given in where, for instance, non-algebraic invariant 
curves with an algebraic cofactor for a polynomial system of degree 4 are presented. 
In our work, we give other families of systems with such invariant curves. 

C. Christopher, in [0], studies the multiplicity of an invariant algebraic curve and 
gives the definition for exponential factor, which is a particular case of invariant for 
system ©• 

Definition 2 Let h, g be two coprime polynomials. The function e h ^ 9 is called an 
exponential factor for system HJ) if for some polynomial k of degree at most d — 1, 
where d is the degree of the system, the following relation is fulfilled: u Ad (e 71 / 9 ) = 
ke h ^ 9 . As before, we say that k(x,y) is the cofactor of the exponential factor e h/ ' 9 . 

We note that an exponential factor e h ^ 9 does not define an invariant curve, but the 
next proposition, proved in j^j, gives the relationship between both notions. 

Proposition 3 J^j If F = e h ^ 9 is an exponential factor and g is not a constant, 
then g = is an invariant algebraic curve, and h satisfies the equation uo A dh = 
(hk g + gk F ) where k g and k F are the cof actors of g and F, respectively. 

All these previous results are closely related to a result of M. F. Singer [23] 
which represents an important progress in the resolution of the integrability problem 
when considering first integrals for a system (J2J) in the class of Liouville functions. 
Roughly speaking, we can define a Liouville function or a function which can be 
expressed by means of quadratures, as a function constructed from rational functions 
using composition, exponentiation, integration, and algebraic functions. A precise 
definition of this class of functions is given in [25] . 

Theorem 4 |25j Let us consider the polynomial 1-form to = Qdx — Pdy related to 
system System $2j) has a Liouville first integral if, and only if, uj has an inverse 
integrating factor of the form V = exp j £, where £ is a closed rational 1-form. 

We notice that the conditions on the function V given in this Theorem can be 
restated as du = £ A to and d£ — 0. 

Taking into account Theorem |3J C. Christopher JU] gives the following result, 
which precises the form of the inverse integrating factor. 

Theorem 5 [1(J| If system has an inverse integrating factor of the form exp j £ 
with d£ = and £ = ^dx + £%dy where £j ; i = 1,2, are rational functions in x and 
y, then there exists an inverse integrating factor of system of the form 

V = exp{D/E}Y[Cl*, 

where D, E and the C{ are polynomials in x and y, and U G C. 



4 



Theorem |S] states that the search for Liouvillian first integrals can be reduced to 
the search of invariant algebraic curves and exponential factors. 

In |2T], H. Poincare stated the following problem concerning the integration of 
an equation ([T}: Give conditions on the polynomials P and Q to recognize when 
there exists a rational first integral. As the same H. Poincare noticed, a sufficient 
condition to solve this problem consists on finding an upper bound for the degree 
of the invariant algebraic curves for a given system (J2J). From Darboux's result, it 
is known that for every polynomial vector field, there exists an upper bound for 
the possible degrees of irreducible invariant algebraic curves. However, it is a hard 
problem to explicitly determine such an upper bound. Some bounds have been 
given under certain conditions on the invariant curves, see D. Cerveau and A. Lins 
Neto work [3], or on the local behavior of critical points, see M. Carnicer's work jl]. 

In this sense, A. Lins Neto conjectured [TH] that a polynomial system (J2J) of 
degree d with an invariant algebraic curve of degree high enough (where this bound 
only depends on d) would have a rational first integral. This conjecture has been 
shown to be false by several counterexamples. J. Moulin-Ollagnier [19 j gives a 
family of quadratic Lotka-Volterra systems, each with an invariant algebraic curve 
of degree 2£, where I is the parameter of the family, without rational first integral. 
A simpler example is given by C. Christopher and J. Llibre in jTIJ. In [S] a family 
of quadratic systems with an invariant algebraic curve of arbitrarily high degree 
without a Darboux first integral nor a Darboux inverse integrating factor is given. 
All these counterexamples exhibit a Liouvillian first integral. 

The natural conjecture at this step, also given by A. Lins Neto, see [T7j . after 
the counterexample of J. Moulin-Ollagnier appeared, is that a polynomial system 
(J2J) of degree d with an invariant algebraic curve of degree high enough (where this 
bound only depends on d) has a Liouvillian first integral. 

In this work we show a relationship between solutions of a class of systems (j2J) 
and linear homogeneous ordinary differential equations of order 2 of the form 

A 2 (x) w"(x) + A x {x) w'(x) + A (x) w(x) = 0, (4) 

where x G R, w'(x) = dw(x)/dx and w"(x) = dw'(x)/dx. We only consider equa- 
tions flU where A { (x) E R[x] for < i < 2 and A 2 {x) 0. 

By means of a change of variable we rewrite an equation (@J) as a polynomial 
differential system such that it has an invariant related to w(x). In case w(x) is a 
polynomial we get an invariant algebraic curve. 

Moreover, we give an explicit first integral for all the systems built up by this 
method by means of two independent solutions of equation (j3J). 

We give analogous results for a linear homogeneous ordinary differential equation 
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of order 1 such as 

w'(x) + A{x)w(x) = 0, (5) 

where x G R, w'(x) = dw{x)/dx and A{x) G K(x). All these results are given in 
Section |21 

In Section El we consider all the families of quadratic systems with an algebraic 
curve of arbitrarily high degree known until the moment of composition of this 
paper and we show that they all belong to the construction explained in Section |21 
The families of quadratic systems with an algebraic curve of arbitrarily high degree 
studied in this paper are the ones appearing in [7J |HJ El EI] and one example more 
first appearing in this work. This new example consists on a biparametrical family 
of quadratic systems, which we give an explicit expression of a first integral for, such 
that when one of the parameters is a natural number, say n, the system exhibits an 
irreducible invariant algebraic curve of degree n. 

We give the explicit expression for the first integral of a certain system (J2J) by 
means of invariant functions for it, and applying the Generalized Darboux's Theory 
as explained in ^3] where a new kind of first integrals, not only the Liouvillian ones 
as in classical theories, is described. We exemplify this result with the families of 
systems depending on parameters described in Section El We remark that the first 
integrals that we give in Section |21 are not, in general, of Liouvillian type. However, 
these first integrals are Liouvillian at the values of parameters which correspond to 
the systems with algebraic solutions. In the Subsection 13.41 we give an example of 
a 3-parameter family of quadratic systems with a center at the origin which can be 
constructed by the method appearing in Section |21 from an equation (JSJ). 

A question suggested by these examples is whether there are polynomial systems 
which are not reversible nor Liouvillian integrable which have a center and can 
be integrated by means of Theorem |H1 see Section 2. The work [2] is related to 
this question as it gives an example of an analytic system, not polynomial, with a 
center which is not reversible nor Liouvillian integrable. All the known families of 
polynomial vector fields with a center at the origin are either Liouvillian integrable 
or reversible, see EH] for the definition of reversibility. In fI7\ |2"H] . Zoladek 
classifies all the reversible cubic systems with a center. The reversible systems may 
have a first integral not given by Liouville functions or no explicit form of a first 
integral may be known. For instance the reversible system x = —y + x 4 , y = x has 
a first integral composed by Airy functions, see |14j . and no Liouvillian first integral 
exists. The system x = —y 3 + x 2 y 2 /2, y = x 3 is an example given by Moussu, 
see [20 , which has a center at the origin since it is a monodromic and reversible 
singular point and no explicit first integral is known for this system. 

Since some examples of polynomial systems, which can be integrated by the 
method described in Section 2, appear after a birrational transformation, another 
suggested open question is if all the polynomial systems with a center are birra- 
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tionally equivalent to one derived from Theorem El or from Theorem ITTl 



2 Homogeneous linear differential equations of 
order < 2 and planar polynomial systems. 

Let us consider a homogeneous linear differential equation of order 2: 

A 2 (x) w"(x) + Ax(x) w'(x) + A (x) w(x) = 0, (6) 

where w'(x) = dw(x)/dx, w"(x) = dw'(x)/dx, Ai(x) G M.[x], i = 0,1, 2, and A 2 (x) ^ 
0. 

Theorem 6 Given g{x,y) = g (x,y)/g 1 (x,y) with gi{x,y) E R[x,y], gi{x,y) ^ 
and dg/dy ^ 0, each nonzero solution w(x) of equation (G|) is related to a finite 
number of solutions y = y(x) of the rational equation 

dy = A (x) gj + A 1 (x) 9l g + A 2 (x) gj + A 2 (x) ( 9l || - g 

dx A ( „ dgi dg \ 

A 2{ x ) ^-3i^J 

by the functional change dw/dx = g(x, y)w(x), which implicitly defines y as a func- 
tion of x. 

Proof. Let us consider equation (jUJ) and the functional change dw/dx = g(x,y)w(x) 
where y = y(x), that is, y is implicitly defined as a function of x. This change may 
also be written as w(x) = exp(f* Q g(s,y(s))ds), where xo is any constant, and it is 
injective. We see that it is not necessarily bijective unless the maximum degree of 
gi(x,y) and g (x,y) in the variable y equals 1. But it defines a finite number of 
functions y(x). 

By this functional change, equation (0) becomes 

w(x) (a q (x) +gA 1 (x) +g 2 A 2 (x) + A 2 (x)^-^- + A 2 (x)^- ) =0. 
\ dx ay ox J 

We have that w{x) is a nonzero solution of © so this equation is equivalent to the 
ordinary differential equation of first order (|7J). Therefore, each non-zero solution 
w(x) of equation (jUJ) corresponds to a finite number of solutions y = y(x) of the 
planar polynomial system ((7j). | 

Theorem 7 We consider the I- form related to equation 

A (x) gf + A^x) g! g + A 2 (x) g% + A 2 (x) (gr^; dx 
„ / n ( dgt dg \ , 



Let w(x) be any nonzero solution of equation (OJ). Then the curve defined by 
f(x,y) = 0, with f{x,y) := gi(x,y)w'(x) — g (x,y)w(x) is invariant for system 
and has the polynomial cof actor 

dgi A ( v dg \ , . dg 

dg 1 dg dg dgi 



k{x,y) 



+A 2 (x) 



dy dx dy dx 



Proof. Let us consider f(x,y) as defined above and let us compute uj A df: 

uj Adf = \^ w \ x ) + 9iw"(x) - ^ w ( x ) - g w'(x) \ A 2 (x) \9q-^ 

dgo\ fdgi ,. . dg A r 2 
~ 9l ~dy J \~dy W ~ ~dy W J ^ 9l + Ai y x )9i9o 



+A 2 {x)gl + A 2 {x) [ gi 



dg dgi 



dx 



go 



dx 



, 9^1 %A , / X „, X . 

g-i 5f + A (x) g x ) (^-w'(x) ~ ~^ W ^J + 

\ ox ay ox ay J ay J 

Since w(x) is a solution of (JSj), we can substitute y4 2 (x)w"(x) by — Ai(a;)w/(a;) — 
A (x)w(x). Therefore, 



uj Adf 



M x )^ + Mx)—j g 1 + 

A (x) — I A (x) ( — — — 991 
dy V <9y <9x 9y cte 



f(x,y). 



Then, we have that the function f(x,y) is an invariant for system J2J) and has the 
written polynomial cofactor. | 

Theorem 8 Let {wi(x) , w 2 (x)} be a set of fundamental solutions of equation (OJ). 
We define fi(x,y) := gi(x, y)w' i (x) — go(x, y)wi(x), i — 1,2. T/ien, system /ias 
a /irst integral H(x, y) defined by 

fifay) gi(x,y)w' 1 (x) - g (x,y)w 1 (x) 



H(x,y) :-- 



h{x,y) g x {x,y)w' 2 {x) - g (x, y) w 2 {x) ' 
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Proof. By Theorem[71 we have that y) := gi(x, y)w' i (x)—g (x, y)wi(x), i = 1,2, 
are invariants for system (J7J), both with the polynomial cofactor 

u \ (a ( \ d 9i , a ( \ d 9o\ ,A f \ d 9o , A ( s. fdgx dg dg dg 
k(x,y) = A (x) — +A 1 {x) — g 1 +A 2 (x) g — — hA 2 (x) 



ay ay J ay \ ay ox ay ox 

We remark that /1//2 cannot be constant since the two solutions Wi(x), i = 1,2, 
are independent. Therefore, 

a j tt f2(uJAdf 1 )-f 1 (coAdf 2 ) f 2 kf l -f 1 kf 2 _ n 
to A ah = - 2 = - 2 = 0. 

J 2 J 2 

So, H(x,y) is a first integral of system (jZJ). | 
Lemma 9 TTie function defined by 

( f x Ai(s) 
q{x) := A 2 (x) exp / -ry^cfe 
VAo M s ) 

is an invariant for system with cofactor (Ai(x) + A' 2 {x)) (go^ ~ di^jf 

We notice that q(x) is a product of invariant algebraic curves and exponential factors 
for system with complex exponents. 
Proof. We compute u A dq and we have 

+A' 2 (x) ( dg x dg Q 



uAdq = w A — qdx = (A x {x) + A' 2 (x)) g - g 1 — q. 

A 2 (x) \ dy dy J 

We notice that this algebraic cofactor has degree < d — 1 provided that system (J7j) 
has degree d. | 

Proposition 10 We use the same notation as in Theorem^ Letw(x) be a nonzero 
solution of and we define f(x, y) := w'(x) — g(x, y)w(x) and q(x) as in Lemma 
The function V(x,y) = q(x)f(x,y) 2 is an inverse integrating factor of system 

Proof. We only need to verify that u A dV + Vduo = 0. We have that 

du= - 2A (x)g 1 — + A 1 (x)\g — + gi—)+2A 2 (x)g -Q- 

+2A 2 (x) ( dQl 990 - dg ° d 9A + A'(x) (g — - gi— ^ 
\ dy dx dy dx J \ dy dy J 
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Then, 

uj A dV 



u A{2qf df + f 2 dq) = 2qf {u Adf) + f 2 {u A dq) 

2 A {x) g 1 — +2Ai{x)g 1 — +2 A 2 (x) g -j- 
dy ay ay 



+2A 2 (x) 



dgidg dg Q dgA , . , N dg x dg 



dy dx dy dx 



+ A 1 (x)g -^--A 1 (x)g 1 ^ 



i At / \ f d 9i dg 
+A 2 (x) \^—- 9l — 

2A (x)g 1 ^ + A 1 (x) 
dy 



V 

dgi dg \ . oA ,s dg 
go^- +9i-z- + 2 M\x) go -Tr- 
oy dy J dy 



+2A 2 {x) 
-Vdw. 



dgi dg dg dg 



dy dx 



+ A'(x) (g — -gi — 
dy dx J \ dy dy 



V 



We remark that Theorem |H1 gives, in general, a non Liouvillian first integral for 
the planar polynomial systems (J7jl . In Section El we analyze some polynomial sys- 
tems constructed from Theorem |H1 that have no Liouvillian first integral. 

We consider now a linear homogeneous ordinary differential equation of order 1 
such as 

w'(x)+A(x)w(x)=0, (8) 

where x E R, w'(x) = dw(x)/dx and A{x) = A (x)/Ai(x) with Ai(x) E M[x] and 
A\(x) 0. We give analogous results for this case whose proofs are not given to 
avoid non useful repetitions. 

Theorem 11 Given g(x,y) = g (x,y)/ gi (x,y) with gi(x,y) E M[x,y], gi(x,y) 
and dg/dy ^ and h(x) = ho{x)/h\{x) with hi(x) E M[x] and hi(x) ^ ; each 
nonzero solution w(x) of equation (QJ) is related to a finite number of solutions 
y = y(x) of the rational equation 



dy _ A 1 {x)h {x)g 2 1 -A (x)h 1 {x)g g l -A 1 (x)h 1 {x)(g 1 ^-g ^) 



dgi 



dx 



9o-§r 



A 1 {x) h(x) [g^ 

by the functional change 

g(x,y) — exp ^— J A(s)ds^j J exp (^J A(r)dr ) h(s)ds 



(9) 



W(X) 
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Theorem 12 We consider the I- form related to equation (OJ) 

(■/./■ 



Ai(x) h (x) g\ - A (x) h x {x) g gi - A x (x) h^x) ( gj^- - g dlJl 



- A 1 (x) h x (x) [g x ^ - g ^j dy. 



OX OX 



Let w(x) be any nonzero solution of equation (0), that is, for Cel - {0} we have 
w(x) = C exp (— J* A(s)ds) . 
Then, the function 

f(x,y):=giw(x)-g + giexp^--J A(s)dsj J exp A(r)drJ h(s)ds 



is 



invariant for the polynomial system (0J) ; with the polynomial cof actor 
k(x, y) = -A Q (x) hi(x) #1 + Ai(x) h (x) g x ^ 



+A 1 (x) hAx) 



dgo_dgi_ _ dgi dg 
dy dx dy dx 



Lemma 13 The function q(x,y) = g%(x, y) exp (Jj^ — A(s)dsj is an invariant for 
system (CJ) with the same polynomial cof actor as f(x,y). 

Theorem 14 We use the same notation as in Theorem EH and Lemma [731 The 
function H(x,y) defined by H(x,y) := f{x,y)/q(x,y) is a first integral for system 
(CJ) and the function V(x,y) := A\(x) h\ (x) g\ (x, y) q(x,y) is an inverse integrating 
factor. 

We remark that H (x, y) is a Liouvillian function and, therefore, a system (jUJ) has 
always a Liouville first integral. 

In Section [3] we give an example of a 3-parameter family of quadratic systems 
with a center at the origin which can be constructed following Theorem ^2 

3 Examples of families of quadratic systems. 

3.1 Quadratic systems with invariant algebraic curves of 
arbitrarily high degree linear in one variable. 

We first consider the examples of families of quadratic systems with algebraic so- 
lutions of arbitrarily high degree appearing in In that work all the invariant 
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algebraic curves linear in the variable y, that is, defined by f(x, y) = pi(x)y+p 2 (x), 
where p\ and p 2 are polynomials, are determined. 

The example appearing in [8j is a further study of an example appearing in 
[7j and the example given in [TT] is also described in [7j. We show that all these 
quadratic systems, with an invariant algebraic curve of arbitrary degree can be 
constructed by the method explained in the previous section. Moreover, we give 
the explicit expression of a first integral for any value of the parameter n, even 
in the case when n is not a natural number. If n is a natural number, we obtain 
the invariant algebraic curves of arbitrary degree and a Liouvillian first integral. 
However, when n G" N we obtain polynomial systems with a non Liouvillian first 
integral. 

As it is shown in [7j, all these families of systems can be written, after an affine 
change of variables if necessary, in the form 

x = n 1 (x), y = (2n + 1) L\x) Q l (x) - ^± ^ Sl^x) Sl'Kx) - L(xf + y\ (10) 

where Qi(x) is any quadratic polynomial, L(x) is any linear polynomial and ' = 
d/dx. We have that system (jTUJl has an invariant curve f(x, y) = 0, where f(x, y) : = 
Pi(x)y + Qi(x)p' 1 (x) — L(x)px(x), with a cofactor y + L(x), where pi(x) is a solution 
of the second order linear differential equation 

T) 

n x {x)w"{x) + {^[{x) - 2L{x))w'{x) + -(4L'(x) - {n + l)n'({x))w{x) = 0. (11) 

In [7] it is shown that, in case n G N, an irreducible polynomial of degree n be- 
longing to a family of orthogonal polynomials is a solution of equation (JUJ. For 
instance, when Q\(x) = 1, we get the Hermite polynomials, when Qi(x) = x, we get 
the Generalized Laguerre polynomials and when Q\(x) = 1 — x 2 , we get the Jacobi 
polynomials. 

We consider again the general case in which n£l and we define A 2 (x) := Qi(x), 
A x (x) := Q[(x) - 2L(x) and A (x) := f (4L'(x) - {n + We have the linear 

differential equation (jllj) in the same notation as in Theorem |H1 and we consider 
L(x) - y 

The system obtained by the method explained in Section El exactly coincides 
with system (fTUj) . We consider a set of fundamental solutions of equation (fTTj) 
{wi(x), W2(x)} and applying Theorem |Hl we have a first integral for system (flT)|) for 
any value of the parameter n G R. 

In case n G N we have that wi{x) degenerates to a polynomial and w[(x) — 
g(x,y)wi(x) = coincides with the algebraic curve given in the work 
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We explicitly give the first integral for each of the families described in [7j and 
for net. We have that A 2 (x) is a non-null quadratic polynomial in this case, and 
depending on its number of roots, we can transform it by a real affine change of 
variable to one of the following forms: A 2 (x) = 1, x, x 2 , 1 — x 2 , 1 + x 2 . 

If A 2 (x) = 1, we can choose L(x) = x by an affine change of coordinates. A set 
of fundamental solutions {wi(x), W2(x)} for (fTTj) with n G K. is 



wi(x) = 2 n v^F 
w 2 {x)= 



r(^) 
r(¥) 



n 1 

'r 2 ; 



n 1 



x 



2x 



r(-|) 

2.T 



1^1 



1^1 



1-n 3 



2 '2' 



1 -n 3 



.r 



2 ; 2 ;3; 



where is the Euler's-Gamma function defined by T(x) = f^°t x 1 e t dt and 
i-F\(a; b; x) is the confluent hypergeometric function defined by the series 



L Fi(a;6;x) = ^ 



k=0 



{a)kX k 
(b) k k\ 



with (a)^ = a(a + l)(a + 2) . . . (a + k — 1), the Pochhammer symbol. See [JJ for 
further information about these functions. So, a first integral for this system is the 
expression given in Theorem |H| H(x,y) := fi(x, y)/f 2 (x, y), where 



f\,2{x,y)= ±r 



1 — n 



• ) I 1 Tb 3 o 

6 (xy - x +1) 1F1 ( — ^— ; -; x 



-4(n- l)x 2 iF x 



3 — n 5 



2 '2' 



+ 



n 

3FI-- 



2nx i-Fi ( 1 - -; -; x 2 J + (x - y) iF x f --; -; x 2 



When n G N, we have that (jllj) corresponds to the equation for Hermite polyno- 
mials and w\{x) coincides with the Hermite polynomial of degree n. The invariant 
algebraic curve given in [7j corresponds to fi(x,y) = 0. 

If A 2 (x) = x, we choose L( ), where a is an arbitrary real constant, 

and a set of fundamental solutions for (fTTj) is: 



Wi[X J 



p/ + j-j l-fi (-w; a + 1; x) , w 2 (x) = x a iFi (-« - n; 1 - a; x) . 
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The first integral for this system is H(x,y) = x a hi(x,y)/h 2 (x,y) with 

h\(x, y) = (2y — x + a) (a + 1) 1F1 (— n; a + 1; x) — 2nx \F 1 (1 — n; a + 2; x) , 

h 2 (x, y) = (2y — x + a) (a — 1) iFi (—a — n; 1 — a; x) 

—2 (a + n) x \F\ (1 — a — n; 2 — a; x) . 

The first integral as given in Theorem |S] is f\(x, y)/f'2{x, y) and we notice that 
H(x,y) = cfi(x,y)/ f 2 (x,y) where c G R — {0}. We do not write c in terms of the 
parameters of the system to simplify notation. 

When n 6 N, we have that (fTTj) corresponds to the equation of Generalized La- 
guerre polynomials and W\(x) coincides with the Generalized Laguerre polynomial 
Ln . The invariant algebraic curve given in [7] corresponds to f\[x : y) = 0, where 
fi(x, y) := w[(x) - g(x, y)w 1 (x). 

If A 2 (x) = x 2 , the birrational transformation yet described in [7j, x = 1/X and 
y = (1/X) (1/2 — Y), makes this case equivalent to the previous one. 

If A 2 (x) = 1 — x 2 , we choose L(x) = + [3)x + (a — /?)), where a, [3 are two 
arbitrary real constants, and a set of fundamental solutions for (fTTj) is: 

(a + l) n / 1 — x 

wi(x) = — -T- 2 Fi — n, 1 + a + [3 + n; a + 1; 



I> + 1) z 1 V 2 

_ ( 1 — x 

w 2 (x) = (1 - a;) - " 2 F X ( -a - n, 1 + (3 + n; 1 - a; — - — 

where 2-F1 (ai, 02; b; x) is the hypergeometric function defined by 

{a 1 ) k (a 2 ) k x k 



2 F X (a u a 2l b;x) = 



k=0 



(&)* k\ 



The first integral given in Theorem |H1 is H(x,y) = (1 — x) a hi(x,y)/h 2 (x,y), 
where 

ht = n (1 + a + (3 + n) (x 2 - 1) 2 Fi ^1 - n, 2 + a + (3 + n; 2 + a; ^-^J + 

/ 1 — x 

(a + 1) ((a + + (a - (3) - 2y) 2 F ± ( -n, 1 + a + (3 + n; 1 + a; 

/i 2 = (a - 1) ((a - + (a + 6) + 2y) 2 F X (^-a - n, 1 + + n; 1 - a; ^—^ ) + 

1 — x 



(a + n) (1 + (3 + n) (x 2 - 1) 2 F X ^1 - a - n, 2 + (3 + n; 2 - a; 

14 



The first integral as given in Theorem |S] is f\(x, y)/f'2{x, y) and we notice that 
H(x,y) = cfi(x, y)/ f 2 (x, y) where c G R — {0}. As before, we do not write c in 
terms of the parameters of the system to simplify notation. 

When n G N, we have that (fTTJl corresponds to the equation of Jacobi poly- 
nomials and w\(x) coincides with the Jacobi polynomial Pn (x) and the invari- 
ant algebraic curve given in [7j corresponds to f\(x,y) = 0, where fi(x,y) : = 
w' x {x) - g(x,y)w 1 (x). 

If A^(x) = 1 + x 2 the complex affine change of variable x = iX makes this case 
equivalent to the previous one, as it is shown in 

We have re-encountered by this method all the examples appearing in [7j from a 
unified point of view. In addition, in this work we have given an explicit expression 
of a first integral for each case and for any value of the parameter n G R. To this end, 
we have found invariants for the system and we have applied the generalization of 
Darboux's method as explained in jT3] to be able to construct a first integral which 
is, in general, of non Liouvillian type. 

3.2 A Lotka-Volterra system 

As it has been explained in the introduction, the first counterexample to Lins Neto 
conjecture was given by J. Moulin- Ollagnier in 19j. His example is a quadratic 
system with two invariant straight lines and an irreducible invariant algebraic curve 
f(x,y) = of degree 11 when I G N. This gives a family of systems depending on 
the parameter i which have a Darboux inverse integrating factor when I G N but 
no rational first integral. The method used in only shows the existence of such 
invariant algebraic curve but no closed formula to compute it is given. We give an 
explicit expression for an invariant by means of Bessel functions for any value of 
t G R — {|} which, in the particular case leN degenerates to the algebraic curve 
encountered in [T^j . 

We show that after a birrational transformation, this example coincides with a 
system constructed with the method explained in Section [21 A birrational transfor- 
mation is a rational change of variables whose inverse is also rational. This kind of 
transformations bring polynomial systems such (J2J) to polynomial systems and do 
not change the Liouvillian or non Liouvillian character of the first integral. 

Let us consider the system appearing in [TH] but assuming that £ G R — {|} 

X = x(l--+y), y = yl-- r - + --yj. (12) 

We make the birrational transformation 

Auv 1 - 2£ 
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whose inverse is 

1 - 2t 

u = xy, v = . 

Ay 

By this transformation, system (|T2|) becomes 

2u 1 - 21 2t + 1 2u 2 

U= i^2t v = —r + 2i^i v + 2i^i v - (13) 

We notice that the equation for the orbits satisfied by the variable v as a function 
of u is a Ricatti equation. 

Let us consider g(u,v) := v and the linear differential equation of order 2 given 

by 

u w"{u) + 1(1 + 21) w\u) -hi- 2tf w(u) = 0. (14) 
2 8 

Applying the method given in the previous section, this linear differential equation 
gives system (JT3j) modulus a change of time. 

A set of two fundamental solutions for equation (JT3jl is given by 

Wl (u) = u^%_ e ((1 - 2£)^£) ,w 2 (u) = u^/ H ((1 - 2l)Jfj , (15) 

provided that i is not of the form |(1 — 2r), with r an integer number, because 
in this case W\ and are linearly dependent. The function I v {u) is the Modified 
Bessel function defined by the solution of the second order differential equation 

u 2 w"(u) + u w'(u) — (u 2 + v 2 ) w(u) = 0, (16) 

and being bounded when u — > in any bounded range of arg (u) with Re(ix) > 0. 
See for further information about this function. 

Hence, by Theorem|H]we have that H(u,v) = fi(u,v)/ f 2 (u,v), where fi{u,v) : = 
w'iiu) —vwiiu) for i = I, 2, is a first integral for system For a sake of simplicity 
we consider the following renaming of the independent variable u = 2z 2 /(l — 2£) 2 . 
This is not a birrational transformation and that's why we only use it to simplify 
notation. The function H writes as: 

(1 - 2£) 2 I (1±2i - ) (z) - 4vzl(2i=i\ (z) 

TT _ \ 2 J \ 2 ) 

{l-2tyi_^^{z) -Avzlfu-,^)' K ' 

By Theorem [7| we have that fi(u,v), % — 1,2 are invariants with the same 
polynomial cofactor k for system (JT3J), so the curve f(u,v) = given by f(u,v) = 
TTZ 2£+1 (f 2 (u,v) — f 2 (u,v)) is also an invariant. We multiply by n only for esthetic 
reasons. 
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Now we assume that i G N and we show that / = defines an invariant algebraic 
curve. To this end we use the following formulas appearing in ^ |2B|. When 
v — \ = n G Z, we define c(n) = —mr^—1/2 and the following relation is satisfied: 



/„(*)= -^e< n \[^\smh(c(n)-z) V 



fc=0 

-I" ' 

L 4 



(M 


+ 2fc- 


\)\ 


(2*)!(|i/| 


-2Jfe- 


§)!(2*)* 



2|r .! \ 



+ cosh (c(n)-z) g (2fc+1)!(M _ 2fc _| )!(2 ,p + i 

where |_^J stands for the greatest integer k such that k < x and \u\ stands for the 
absolute value. 

From the former equation we obtain the following two equalities, with v — | = 
nGZ and leN, 

, 2 

- /_ (/ _|)(«) = (-1) — 

L|J //. . / \ 2i\ A -J 2 J //, -\i / 1 \ 2j'-l 



v (l + 2»)! * (£ + 2j)\ 



_ o (2i)!(^-2i)« V2z7 I \j^o (2j + l)!(^-2j-2)! ^ 



L - ^ //i . ^- . -.m / -, \ 2i-l\ AM" J 



(£ + 2z + i)! f J_Y /n 

^ (2i + !)!(£ - 2i - 1)! \2z) J I ^ (2j)!(£ - 2j - 1)! ^ 



(20) 

Then, we have that fi(z,v) = (1 — 2£) 2 I i+ i(z) — 4vz I £ _i(z) and f2(z,v) = (1 — 
2^) 2 1_( e+ i)(z) — Avz I_( £ _i)(z), and we write / arranged in powers of v: 



f(z,v) = 7 r^ +1 ((l-2£) 4 (/ £ 2 + ,(^)-/! ( , +|) (,)) 

- 8vz(l - 2£) 2 (I iH (z)I e _^z) - I_ {e+ i 2) (z)I_ {e _ L2) (z)) 

Let us consider each coefficient of v in f(z,v) separately and we will show that it 
is an even polynomial in the variable z. The coefficient in f(z,v) of v is: 

vr^ +1 (l-2£) 4 (/ 2 +| (^)-/ 2 (£+|) (z)), 
17 



which by equation (|19|) is an even polynomial in the variable z of degree 21. The 
coefficient in f(z,v) of v 2 is: 



which also by equation ()19|) is an even polynomial in the variable z of degree 2£ + 2. 
Finally, the coefficient in f(z,v) of v 1 is: 

8tt(1 - 2£) 2 z 2i+2 (I e+ ,(z)I e _,(z) - J_ (l+ i ) (z)J_ (l _i ) (z)) ) (21) 

which by equation (}2T)|) is an even polynomial in the variable z of degree 21. 

Hence, we have that f(z, v) is an even polynomial in the variable z of total degree 
21 + 4. When rewriting z = (1 — 2£)y/u/ y2 we have that f{u,v) is a polynomial 
of total degree £ + 2 which is irreducible. The fact of being irreducible is easily 
seen because it is a polynomial of degree two in v and it cannot be decomposed in 
linear factors (the discriminant is not a polynomial raised to an even power) and 
the coefficients of v and v 2 do not have any root in common. 

Undoing the birrational change of variables we get that f(x,y) is an irreducible 
polynomial of degree 2£ given by: 



f(x,y) = x i+1 *y £ -* [2y (lf H (z) - I 2 _ {e+1 Jz)) + 
2V2^ (l i+ i(z)I e _i(z) - I^+l 2) (z)L 
x (ll k {z) - H^{z)) 



where z is the same variable as before, that is, z = — ^ v ^ - 

By equation (fTTjl we can write the first integral for system (fT2"j) for any value of 
£ E R- {±(1 -2r) | r E N}: 



\/2y I(l±2t)(z) - \fxl( 21-11 

H(x,y) = 



2yl /i+2t\(z) - ^Jxl /2t-x\(z) ' 

We have studied system (|T2"j) for any value of the parameter £ e R — {|(1 — 2r) | 
r G N} giving an explicit expression for a first integral using Theorem El and the 
Generalized Darboux's theory as explained in This first integral is not of 

Liouvillian type. Moreover, we give one of its invariants with a polynomial cofactor. 
In the particular case £ £ N, this invariant is the invariant algebraic curve whose 
existence was proved in |19j . 
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3.3 A new example of a family of quadratic systems with 
an invariant algebraic curve of arbitrarily high degree 

We give another example of a family of quadratic systems with an irreducible invari- 
ant algebraic curve of degree 2£ when £ 6 N, where £ is a parameter of the family. 
This family also depends on the parameter a e R. 
Let us consider the quadratic system 

x = {2a - l)£x - a{2£ - l)y + 2a(a - £)(2£ - l)x 2 - 2a 2 (2£ - l) 2 xy, 

y = y(2(2a-l)£+2a(2a-2£-l)(2£-l)x-Aa 2 (2£-l) 2 y), ^ 

where a,!6l which satisfy a ^ 0, £ ^ \ and {2£— l)a 2 — 2£ ^ 0. An straightforward 
computation shows that system (|2*2*Jl has y — and y — x 2 = as invariant algebraic 
curves. 

Let us consider the following birrational transformation x = Y, y = XY 2 whose 
inverse is X = y/x 2 and Y = x. In these new variables system (|22|) becomes 

X = 2a(2£- 1)(X-1)XY. 

(23) 

Y = ((2a- l)£ + a{2£- \){2a ~ 2£ - X)Y - 2a 2 {2£ - l) 2 XY 2 )Y. 

By a change of the time variable we can divide this system by Y and the resulting 
system coincides with the one described in Theorem |B] taking A 2 {X) := 2X(X — l) 2 , 
A X {X) := (2£-2a+3X)(X-l), A (X) := £(l-2a) and g(X,Y) := a(2£-l)Y/(X - 
1). The equation A 2 (X)w"(X) + A 1 (X)w'{X) + A (X)w{X) = has the following 
set of fundamental solutions in this case: 

w x {X) = (X -l)~ e 2 F 1 (l-£,-£;a-£;X 



2 



3 



w 2 {X) = {X -ly^X 1 -^ 2 F X U- a, ^-£-2-a + £-X 

By Theorem [7[ fi(X, Y) = w'^X) — g(X, Y)iVi(X), i = 1,2, define invariants with 
a polynomial cofactor for system (}2"3"|) . Moreover, by Theorem |H] we have a non 
Liouvillian first integral given by H(X, Y) = f^X, Y)/f 2 (X, Y). 

In the particular case £ G N, we notice that f\ (X, Y ) = is a rational function. 
It is an easy computation to show that this rational function is a polynomial when 
rewritten in coordinates x and y. This polynomial gives place to an invariant 
algebraic curve of degree 2£ for system (}2*2*|) . That is, by undoing the birrational 
transformation, we deduce that j\ (x, y) is an irreducible invariant algebraic curve 
for system (|22j) . given by: 



h(x,y) = 2(o -£){£+ (2£ - l)ax)x 2 ^ 2 FA l - ' ' 



2 1 ' x 



2 
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+£{2£ - l)x 2e - 3 (x 2 -y) 2Fl (^-e,l-£;l + a-£;^ 

It is easy to see that the polynomial f\ (x, y) has degree 21 and the cofactor associated 
to the invariant algebraic curve fi(x,y) = is 021 — l)((2a — 1) + 4a(a — £)x — 
A(2£-l)a 2 y). 

The first integral for (|2~2~)) is given by H(x,y) = y a ~ e fi(x,y)/h(x,y), where 
h(x, y)= 2 (a - £ - 2) [(a - £ - 1) x 2 + 

(1 - a - a(2£ - l)x) y] x 7 ~ 2a 2 F X (l - a, | - a; 2 - a - £; ^ j + 

(a - 1) (2a - 3) x 5 " 2a (x 2 - y) y 2 F 1 U - a, | - a; 3 - a + V 

We notice that when both a and £ belong to the set of natural numbers, we have 
that h(x, y) = is an invariant algebraic curve different from j\(x, y) = 0. Then we 
have a quadratic system with a rational first integral H(x, y) with arbitrary degree. 

3.4 A complete family of quadratic systems with a center 
at the origin 

In this subsection we give an example of a 3-parameter family of quadratic systems 
with a center at the origin which can be constructed using Theorem ITT1 The family 
encountered corresponds to the reversible case, see [25] , 

The family of quadratic systems depends on twelve parameters, but up to 
affine transformations and positive time rescaling, we get a family of five essen- 
tial parameters. We have taken a system and we have chosen g(x,y) := y 2 , 
h(x) := 2x(d x — 1)/(1 + ax) and A(x) := 26/ (1 + ax), where a, 6, d are real param- 
eters. Using Theorem ^2 we have encountered the 3-parameter family of quadratic 
systems next described. We remark that in spite of the simplicity of the chosen 
polynomials g(x,y), A(x) and h(x), we amazingly obtain the complete family of 
quadratic systems with a reversible center at the origin. We notice that other 
choices of the functions g(x,y), A(x) and h(x) would give place to other families of 
polynomial systems. 

Let us recall that a center is an isolated singular point of an equation (j2J) with 
a neighborhood foliated of periodic orbits. When the linear approximation of an 
equation (J2J) near a singular point has non-null purely imaginary eigenvalues, the 
point can be a center or a focus. To distinguish between these two possibilities is 
the so-called center problem. H. Poincare gave a method to solve it by defining 
a numerable set of values, called Liapunov-Poincare constants, which are all zero 
when the singular point is a center and at least one of them is not null when it is 
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a focus. When these constants are computed from a family of systems, they are 
polynomials on the coefficients of the family. Hilbert's Nullstellensatz ensures that 
there always exists a finite number of independent polynomials which generates the 
whole ideal made up with all these Liapunov-Poincare polynomials. The zero-set of 
these independent polynomials gives place to the center subfamilies. The reader is 
referred to (2211211 f° r a survey on this subject. 

The computation of these center cases for the family of quadratic systems was 
done by Dulac [|13j for the case of complex systems and a proof for real systems 
is given in [T%] . We also refer Bautin [3| who showed the existence of only three 
independent constants. The computation of the zero set of these three independent 
values gives place to four complete families of quadratic systems with a center at 
the origin which are described in [24j. 

Let us now consider an equation @ such as (1 + ax) w'(x) + 2bw(x) = and 
g(x,y) and h(x) as formerly defined. The rational equation as constructed in The- 
orem ^2 is 

dy —x + dx 2 — by 2 
dx y + axy 

which gives the corresponding quadratic planar system 

x = y + axy, y = —x + dx 2 — by 2 . (24) 

We suppose that ab(a + b)(a + 26) (a + b + d) ^ 0. In case this value is zero, the 
origin of system (f2"4*j) is still a center but with a Darboux integrating factor instead 
of a Darboux first integral. This particular case can also be studied by our method, 
but we do not write it to avoid giving examples without essential differences. 

By Theorem 1121 we have that f(x, y) is an invariant of system (j24j) with cofactor 
—2by, where f(x,y) is given by 

f{x,y) = b (a + b) (a + 2b) w(x) — (a + b + d) (1 + ax)~~ 

-b (a + b)(a + 2b) y 2 + b(a + 2b) dx 2 (25) 
-2b (a + b + d) x + a + b + d , 

with w(x) any non-zero solution of (1 + ax)w'(x) + 2bw(x) = 0, that is, w(x) = 
C(l +ax)~^. 

Choosing C = (a + b + d) (b(a + b)(a + 2b))~ 1 we get an invariant conic. System 
(I24J) has two invariant algebraic curves, the former conic with cofactor —2by and 
an invariant straight line given by 1 + ax = with cofactor y. The Darboux first 
integral 

H(x,y) = (1 + ax)~f(x,y) 
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coincides with the first integral described in Theorem EI 

The origin of this system is a center since it is a monodromic singular point with 
a continuous first integral defined in a neighborhood of it. This example addresses 
to the thought that other families of polynomial systems of higher degree with a 
center at the origin can be easily obtained by this method, avoiding the cumbersome 
computation of Poincare-Liapunov constants. 
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